Creation and destruction of a spin gap in weakly coupled quarter-filled ladders 
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We investigate weakly coupled quarter- filled ladders with model parameters relevant for NaV2 05 
using density-matrix renormalization group calculations on an extended Hubbard model coupled 
to the lattice. NaV205 exhibits super-antiferroelectric charge order with a zigzag pattern on each 
ladder. We show that this order causes a spin dimerization along the ladder and is accompanied by 
a spin gap of the same magnitude as that observed experimentally. The spin gap is destroyed again 
at large charge order due to a restructuring of the spins. An analysis of an effective spin model 
predicts a re-creation of the gap by inter-ladder singlets when the charge order increases further. 
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The discovery of a phase transition at Tq « 34if 
in NaV2 05, below which charge order and a spin 
gap appear 0, has precipitated intensive theoretical in- 
vestigation. NaV205 consists of well-separated planes 
which contain weakly coupled quarter-filled vanadium 
ladders 3] (Fig. They can be described by an ex- 
tended Hubbard model (EHM) 0. A zigzag charge or- 
der as observed in NaV205 is already created in this 
model by the nearest-neighbor Coulomb repulsion V for 
an isolated ladder, but only at overly large values_of V 
|ili,0llliE3- In a recent DMRG study we 
showed that the inclusion of a strong effective Holstein 
couplin g to the lattice, which was found in LDA calcu- 
lations [ij, reduces the required Coulomb repulsion to 
a realistic value. The DMRG calculations then yielded 
good agreement between theoretical and experimental re- 
sults for the amount of charge order, the extent of lattice 
distortion in the c-direction, the effective spin coupling 
Joff in the 6-direction, and the charge gap. 

However, on an isolated ladder the spin gap van- 
ishes in the thermodynamic limit. The occurrence of the 
spin gap in NaV205 appears to be intimately connected 
to the coupling of ladders, as indicated, e.g., by the split- 
ting of magnon branches 0, IT3. . An intriguing sce- 
nario has been put forward by Mostovoy and Khomskii 
[l^. It is based on the fact that the experimentally 
observed unit cell of NaV205 in the ordered phase is 
2a X 26 X 4c The charge order then has a period- 

icity of four parallel ladders in the a-direction, and is 
thus super-antiferroelectric (SAP) The correspond- 
ing polarization of the electrons on the rungs is illustrated 
by wide-narrow ellipses in Fig. ^b) (top). Each ladder 
exhibits antiferroelectric charge order along the ladder 
(generally called zigzag charge order) as well as an anti- 
ferroelectric order to the next nearest ladder. The lad- 
ders correspond to effective antiferromagnetic spin- 1/2 
chains in the 6-direction. For each ladder, the two clos- 
est sites on the neighboring in-plane ladders alternate be- 
tween low charge occupation [indicated by dashed loops 
in Fig. n^b)] and large occupation. Large charge occupa- 
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FIG. 1: Quarter-filled coupled ladders in NaV205. The map- 
ping of a rung (top) onto an effective spin site (bottom) is 
illustrated for (a) a charge-disordered and (b) a SAF-ordered 
regime. The shifted periodic boundary conditions in the a- 
direction are indicated by numbers which identify identical 
sites. The proposed singlet formation for NaV205 is depicted 
by dashed ellipses in the effective spin model (b, bottom). 



tion should effectively weaken the electron hopping along 
the ladder, both through purely electronic interactions as 
well as by pushing away the neighboring oxygen atoms 
and thus reduce the spin coupling. With SAF charge 
order, the effective spin chains are therefore dimerized, 
which was proposed to lead to the formation of spin sin- 
glets [Fig. n^b), bottom] and to the observed spin gap. 
This scenario is difhcult to evaluate quantitatively. It has 
been investigated using exact diagonalization on 
small clusters and by means of an approximate xy-model 
together with a mean- field approximation . It has also 
been argued that the mechanism proposed by Mostovoy 
and Khomskii would not work for NaV205 [l8j. 

In the present paper, we use the Density Matrix Renor- 
malization Group (DMRG) to study the extended Hub- 
bard model with coupling to the lattice ^ for large sys- 
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tems of coupled ladders up to length 20, with periodicity 
of four ladders in the a-direction. We show that the sce- 
nario proposed by Mostovoy and Khomskii indeed works. 
In the thermodynamic limit, it produces a dimerization 
and a spin gap of similar size to that seen in the exper- 
iment. The spin gap is found only in the case of the 
four-ladder SAF order; it does not appear for charge or- 
dering with one- or two-ladder periodicity as studied in 
previous papers Surprisingly, the spin gap closes 

again at large Coulomb repulsion V. We explain this be- 
havior by considering an effective spin model for which 
DMRG calculations on larger systems are possible. 

Model. We use the extended Hubbard model with cou- 
pling to the lattice as introduced in Ref. 10, H = i/EHM + 

Hi + Hc-i, where iJEHM = - + ^^-c-) + 

C/X^i^-jT^U + Vij^inj , with hopping matrix ele- 

ments from first-principles calculations 0, 0| , ta = 
0.35ey, tb = 0.5ta, tah = 0.17ta and 0.33ia, and a 
uniform on-site repulsion U = S.Ota- The lattice and 
the corresponding model parameters are illustrated in 
Fig. n^a). Since the nearest-neighbor Coulomb repul- 
sion V is difficult to estimate, we investigate the func- 
tional dependence on this parameter. Recent LDA cal- 
culations have shown that there is a very large 
Holstein-like coupling to out-of-plane movements of oxy- 
gen atoms, which can be modeled by including the terms 
Hi = and H^-i = ZiUi, with k = 0.125 ^a 

and C = 0.35 ta, and Zi in units of 0.05 A, in the Hamil- 
tonian. We treat these movements adiabatically like in 
Refs. ^3 and 0. The optimal configuration for the dis- 
tortions Zi is then a zigzag-pattern 8] with amplitude 
Zopt, for which we take the results from an isolated lad- 
der '20"|. Then, Hi becomes an irrelevant constant and 
Hc-i reduces to a zigzag alternation of the local chemical 
potentials. 

In order to investigate the proposal by Mostovoy and 
Khomskii, a system with a periodicity of four ladders 
is necessary. Within the DMRG, calculations on an ex- 
tended Hubbard model of four coupled ladders would be, 
however, too restricted in length. We therefore employ 
a system of only two coupled ladders of length up to 20, 
but with shifted periodic boundary conditions, as illus- 
trated in Fig. ^ They ensure the proper periodicity and 
the same SAF structure of neighboring charges in the or- 
dered phase as in a four-ladder system. We apply open 
boundary conditions (obc) in the 6-direction. 

Results. We calculate the spin gap using 0, 0| 
As(i) = Eq{L,N,S, = l)-Eo{L,N,S, = 0) , 
where Eq is the ground state energy, L is the length 
of each ladder, N ~ L/2 is the number of elec- 
trons, and is the total spin in z-direction. The re- 
sults are extrapolated using linear and quadratic fits 
in l/L, including L = (8), 12, 14, 16, 20, as illustrated 
in Fig. Ola) . Due to broken translational invariance in 
the 6-direction {obc), we can define a spin dimerization 





0.15 
0.1 
0.05 




^ V=l.3,m^y=0.57, t^^=0.33l 

□ V=l.3, mj,^j=0.57, t^^=O.I7l 

O V=O.9,m„=0,l_j=0.17t_ 

quadratic fit (L=8-20) 

--. linear fit (L=12-20) 




Hubbard model 




Hubbard model 



m„=0.0, jVi|I=1«5 
mj,y=0.9, jVJh=l/45 
quadratic fit (L > 20) 
quadratic fit (L=8-12) 



(b) 



O- 0.4 



-.aP^l^ ■ ■ ■ ■ ■ 

^i?r\^^'^ eff. spin model 







r (d) 



eff. spin model 



- 0.15 

□ : 

- 0.05 








0.05 0.1 
l/L 



0.05 0. 1 

l/L 



FIG. 2: Examples of finite-size extrapolations for (a,b) the 
spin gap As and (c,d) the dimerization d of (a,c) the Hub- 
bard model and of (b,d) the eff'ective spin model. Linear and 
quadratic fits in ^ are illustrated for charge-disordered states 
{mco = 0, circles) and SAF-ordered states (mco > 0). 



d a.s d — 2 4 J2i=i,nJ2ii-^y {^1,1(^1+1,1 ~ ^1-1,1) 
where i counts the rungs in the ^-direction. We re- 
strict to four rungs in the middle of each ladder 
{I — I, II). Here S^^ is the sum of the S'^-spin at the 
i-th rung. The dimerization d is an indicator for bond 
alternation and spin singlets along the ladder, as illus- 
trated in Fig. n^b) (bottom). Translational invariance 
in the 5-direction is additionally broken by the lattice 
distortions, and we then find the contributions to d to 
be positive on both ladders. A state consisting of con- 
secutive singlets in the a-direction [Fig. ^b), bottom] 
would give the maximum value for d, i.e., Idmaxl = 
1/4. The charge order parameter mco is given by [loj 
m 
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Q = (tt, tt) and N is the total number of sites on the 



ladder. We calculate 1 
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for each ladder separately and 



then average. In the thermodynamic limit, we obtain 
the results shown in Fig.|21for the charge order parame- 
ter mcO) the spin gap As, and the dimerization d, as a 
function of V . For comparison, we also show the lattice 
distortion Zopt and the effective spin coupling Xg, which 
were determined for an isolated ladder in Ref. |8l We can 
identify three different regimes: 

(i) For small nearest- neighbor Coulomb repulsion 
V < 0.95 ta, the lattice distortion and the charge order 
are both zero. This regime is therefore not relevant for 
NaV205. The extrapolated values of the spin gap As 
are mostly finite, but vary strongly between linear and 
quadratic extrapolations, reflecting the difficulty of a re- 
liable estimate in this regime. However, in the 2d ther- 
modynamic limit we can expect a vanishing gap as long 
as mco — 0, from calculations on a system with pure pbc 
in the a-direction |^ and from considerations of an effec- 
tive spin model (see below). The dimerization d, which 
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FIG. 3: Results for the EHM with couphng to the lattice, as 
a function of V. Extrapolated values of (a) the charge order 
parameter mco, (b) the spin gap As, and (c) the dimeriza- 
tion d are illustrated for two coupled ladders (tah = 0.17 ta)- 
The effective magnetic exchange constant Jeff [g| and the dis- 
tortion Zopt of an isolated ladder are also shown in (a). 
The three regimes discussed in the text are highlighted. 

is expected to be zero in this regime, shows slightly neg- 
ative values, probably due to finite-size effects. 
(ii) At a critical Vc ~ 0.95 0j the lattice distortions 
become finite and cause charge ordering in a SAF pat- 
tern for our choice of boundary conditions. Linear and 
quadratic extrapolations to L = oo match quite well and 
we find that the spin gap As and the intra-ladder dimer- 
ization d are finite. This is in contrast to the non-SAF 
charge-ordered system of two coupled ladders with purely 
periodic boundary conditions studied in Ref. Hi which 
shows neither a finite spin gap nor finite dimerization 
in the thermodynamic limit. The concurrent appearance 
of intra-ladder dimerization and a spin gap only when 
SAF charge ordering occurs strongly indicates that here 
charge ordering indeed causes singlet formation along the 
ladder. 

(in) At larger V « 1.6 1,, the spin gap vanishes again, 
to our initial surprise [22|. Concurrently, the dimeriza- 
tion becomes smaller, but does not completely disappear. 
This behavior and the negative extrapolations for the 
spin gap are consistent with the effective spin system 
(see below). 

Regime (ii) can be associated with NaV205 because it 
provides a description for the concurrent appearance of 
lattice distortions, zigzag charge order, and the open- 
ing of a spin gap in the low temperature phase. At 
the effective V = 1.3 ta previously determined lat- 
tice distortions, charge ordei\ as well as the spin gap 
(Ag'^P « lOmeV = 0.029 ta 11) are of the same magni- 
tude as in NaV205 Calculations for different inter- 



ladder hopping strengths tab show that the maximal size 
of the gap increases with tab {tah = 0/0.17/0.33 ta 
As, max ~ 0/0.025/0.05 ta), but that concurrently the 
transition from regime (ii) to (Hi) appears at smaller V 
(for tab = 0.33 ta at y « 1.4 ta), corresponding to smaller 
mco- 

Effective spin model. In order to better understand 
the behavior of the spin gap in the EHM, we consider an 
effectiv e sp in model for NaV205 which was introduced 
in Ref. [13. This model allows us to investigate much 
larger systems - up to length L = 120 - and thus to ex- 
trapolate results to the thermodynamic limit more accu- 
rately. By replacing the quarter-filled rungs (one electron 
on each rung) by a single effective spin site, we obtain 
the Heisenberg spin model whose lattice structure is de- 
picted in Fig. da) (bottom) with two different magnetic 
exchange constants, J° along the effective chains, and 
J°^j, inter-chain. The chains in this effective spin model 
correspond to the ladders of the original system. In the 
SAF charge-ordered state with weakly coupled effective 
chains [Fig. ^b)], the magnetic exchange J^^ differen- 
tiates to three different values depending on the posi- 
tions of the electrons on the interacting rungs ■ large 
■J' = '^ab (1 + mco)^, medium J" = J^^ (1 - mco); ^^'^ 
small J'" = jo^ (1 - mco)^ as illustrated in Fig. [Jb). 
The bond alternation along the 6-direction is taken to be 
Ji/2 = •^b (1 ± 5), where Jb ~ J^^ {I — m^o): and we set 
(5(mco = 0) = and 6{mco « 0.44) « 0.034 13]. As- 
suming a linear dependence yields i5(toco) ~ 0.076 mco- 
Together with jojjo = 1/45 ^ (or J°J4 = 1/25), 
where J°,^ and are the exchange constants at mco — 0, 
the parameters of the effective spin model are completely 
determined and the behavior can be investigated as a 
function of the charge order parameter, mco- 

Like in the Hubbard system, we apply the DMRG to 
two coupled spin chains (corresponding to two Hubbard 
ladders) with shifted pbc in the a-direction, and obc as 
well as pbc in the 5-direction. Due to the small inter- 
chain couplings, the finite-size effects in the a-direction 
are small |24| . justifying calculations on only two coupled 
effective chains. A quadratic extrapolation [Fig.[2Ib),(d)] 
of the finite system sizes (L = 20—120 for obc, L = 20—60 
for pbc) to the thermodynamic limit provides the mco- 
dependence of dimerization d and spin gap As. The 
dimerization and the spin gap increase with mco up to a 
maximum (Fig. ; then, both decrease even though the 
bond alternation still increases. This behavior matches 
well with the results from the EHM. In contrast, un- 
coupled chains (ID spin-l/2-Heisenberg chains) show a 
continuing increase of dimerization and spin gap. 

The destruction of the spin gap is likely due to 
a restructuring of the singlets, which now occur with a 
larger probability on the J'-bonds. For small charge or- 
dering, all inter-chain exchanges are of comparable size 
and are small in units of Jb • The singlet formation along 
the 6-direction is then only marginally influenced. On 
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FIG. 4: Effective spin model. Quadratically extrapolated val- 
ues (L = oo) for (a) the spin gap As and (b) the dimerization 
d, as a function of the charge order parameter mco- Cal- 
culations with obc (black dots) and pbc (open circles) in the 
6-direction are shown. For comparison the spin gap As for 
the EHM is shown, in units of Jeff [Fig. I^Ja)]. 



Mostovoy and Khomskii '14'|. This regime provides the 
possibihty to understand the concurrent appearance of 
charge order, lattice distortions, and spin gap in NaV205 
and is in good agreement with experimental observations 
for NaVaOg. 

Higher charge order, unrealistically large for NaVaOg, 
results in an increased inter-ladder magnetic exchange 
and a restructuring of the spins. The correspond- 
ing destruction and re-creation of the spin gap illus- 
trates an interesting transition between different spin gap 
regimes. While regime (ii) is well described by weakly 
coupled rectangular ladders with dimerization inducing 
SAF charge order, regime (Hi) displays weakly coupled 
zigzag ladders with a charge bond density wave |2R| |. 

This work was supported by the Austrian science fund 
FWF, project P15520. We thank M. Aichhorn, C. Gros, 
T.C. Lang, and F. Michel for useful conversations. 



the other hand, a large value of mco drastically changes 
the inter-chain exchanges, e.g., for mco = 0.9 and 
jo^/jo = 1/45, they become J'/Jb « 0.4, J"/ Jb « 0.02, 
and J'"/Jh « 0.001. The large inter-chain exchange 
J' destroys the spin gap and causes a restructuring of 
the singlets. For pbc along the ladder [open circles in 
Fig. 01a)] a re-creation of the spin gap is observed af- 
ter passing a critical charge ordering m^Q. In contrast, 
o&c-calculations (black dots) show a vanishing gap for all 
"mco ^ '^CO' although both calculations match excel- 
lently for mco < "^co [Fig-El^)]. This strongly suggests 
a singlet formation on the J'-bonds Q with increasing 
J'. In the case of obc, unpaired spins are then present 
at the ends of the ladders [Fig. ^b)], preventing a spin 
gap. This explanation is supported by a considerable in- 
crease of inter-chain spin correlations on the J'-bonds for 
TOGO > mco- 

For obc, after passing the critical charge order m^oj 
the spin gap extrapolates to values that are slightly neg- 
ative [Fig. Ela)]. If the extrapolation is restricted to 
smaller system sizes [dotted line in Fig.|2Ib)], the values 
become more negative. This is in agreement with the 
negative extrapolations for Ag in the Hubbard model in 
the large mco-regime [Fig.|2Ib)] . The magnitude of these 
negative values gives a minimal estimate of the error in 
the extrapolation. 

Conclusions. We have studied coupled quarter-filled 
ladders in an extended Hubbard model with lattice 
coupling and with model parameters appropriate for 
NaVaOs. Our calculations match the effective spin model 
proposed by Gros and Valenti [l^ . and show two ways 
to realize a spin gap for coupled quarter-filled ladder 
compounds. At lower charge order we find a super- 
antiferroelectric regime with spin dimerization along the 
ladder and a spin gap, in agreement with the scenario of 
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